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Class of axially symmetric solutions of vacuum Einstein field equations including the 
Papapetrou solution as particular case has been found. It has been shown that the de- 
, rived solution describes the external axial symmetric gravitational field of the source 

» I ■ with nonvanishing mass. The general solution is obtained for this class of functions. 

' As an example of physical application, the spacetime metric outside a line gravitomag- 

netic monopole has been obtained from Papapetrou solution of vacuum equations of 
gravitational field. 
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1. Introduction 

A wide range of massive objects from the compact relativistic stars, as black holes 
and neutron stars, up to modern cosmological models, including Early Universe 
and its evolution 1 where the solutions of Einstein equations are applicable and 
important is currently known. But the logic of scientific development and practical 
experience requires an introduction into the consideration of principally new ob- 
jects, like cosmic strings, gravitomagnetic monopoles, quasiclosed worlds etc. The 
extension of objects under research leads to the additional analysis and extension 
of well known existing solutions of fundamental field equations. Therefore reconsid- 
eration of old and finding new exact exterior solutions of Einstein field equations 
describing space-time around these objects is one of the interesting and important 
branches of general relativity and theoretical physics. There exist a lot of inter- 
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csting problems, especially in quantum gravity, because it is still not clear how to 
build quantum gravity (see, for example, 2 ) and in this level of our understanding 
it is very useful to investigate objects which are on the boundary between classical 
and quantum gravity. Of special importance also are external gravitational fields 
of isolated stars, as today is recognized they are at the core of most intriguing as- 
trophysical phenomena. In spite of the fact that many exact solutions of Einstein 
equations can be found in classical monographs as 3 it is still interesting to study 
new possible solutions of this fundamental equation. 

In our present work we concentrate our efforts in investigation and analysis of 
some exact solutions of Einstein equations 3 , namely, on the exterior axisymmetric 
fields in the vacuum, arising around rotating gravitational objects. In our consider- 
ation we look for some mathematical features and structure of Einstein equations, 
which could determine their physically acceptable solutions. 

In fact, our motivation is, that the Papapetrou class of axially-symmetric so- 
lutions 4,5 of vacuum Einstein equations does not contain physically interesting 
one, because in the classical Newtonian limit there is no term being proportional 
to r _1 , what means that it describes rotation field for the bounded systems with 
zero mass 5 . However it should be noted, that one of the solutions of Einstein equa- 
tions with the similar asymptotic behavior has been considered and interpreted 
as gravitational field of massless quadrupolc 6 . From geometrical point of view, 
the gravitational source is characterized by a topological closure of internal three 
dimensional space with zero integral contribution according to its zero mass. 

In the Papapetrou case the solutions are given in the class of pure harmonic 
functions £ and our assumption is to investigate this solution in the more extended 
functional space of functions £. 

Here we follow to the logic of investigation, brilliantly reflected in 5 . As we will 
show in the Section 2, small difference, introduced by us into the approach described 
in 5 , gives us a possibility to obtain solution of Einstein equations for more general 
class, which includes Papapetrou one as particular case. 

Throughout, we use a space-like signature (+, — , — , — ) and a system of units 
in which G = 1 = c. Greek indices are taken to run from to 3 and Latin indices 
from 1 to 3; partial derivatives are denoted with a comma. 

2. Extension of Papapetrou class of solutions 

We consider a stationary axially symmetric object in cylindrical coordinates x a = 
(t, p, <f), z) and assume that there is no any gravitating matter outside of it. For the 
rotating fields in GR, the most general form of the metric of space-time of axially 
symmetric objects can be written as Weyl - Lewis - Papapetrou one, presented in 
the form 5 

ds 2 = f(dt - ujd(f>) 2 - p 2 • / _1 # 2 - e»{dp 2 + dz 2 ) , (1) 

where fl + k 2 = p 2 , f, u), p, k, I are unknown metric functions may depend on the 
cylindrical coordinates p and z. 
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The external gravitational field is described by the symmetric Ricci tensor R a p 
which obeys to the exterior Einstein vacuum equations R a p = 0. The nonvanishing 
Rqo and i?03 components of Einstein equations can be written as 5 

f(l PP + f,zz + P -\f, P ) f 2 p fl + P~ 2 f\u 2 p + w?J = , (2) 

f(uj. pp + u tZZ - p _1 w, P ) + 2f, P u, P + 2f,zU,z = . (3) 

Following to the approach presented in 5 one can select the function uu satisfying 
to equation 

w,pp + ^.zz - P^u.p = , (4) 

as u = Cp(, p , where C is a constant. Mathematical properties of function ( will be 
defined below. It is easy to see, that after substituting u> in (4) and using simple 
relations 

w p = C(C,p + pC,pp), v,pp = C(2(. pp + pCppp) , 

u,z = CpC.pz, U.zz = Cp(, pzz , (5) 

one can write 

C PPP + C P zz + (p- 1 C P ), P = o . (6) 

Integration of the equation (6) gives an expression for ( : 

V 2 C = C, PP + Lzz+ P - 1 L P = r,(z) . (7) 

Therefore we introduced an arbitrary function rj = r)(z). Introduction of this 
function is very important, because it could generate different kind of solutions for 
set of differential equations (2) and (3). Hereafter we have essential difference with 
compare to 5 and our aim is to solve Einstein equations (2) and (3) in terms of 
more wide class of functions (, which include harmonic functions as a particular 
case. 

In the case r](z) = equation (7) becomes Laplace one, which is similar to the 
expression given in 5 . However in the general case when rj(z) ^ 0, we have Poisson 
like equation (7) for function £ instead of Laplace one. Using the relations (5) for 
the function ui and condition (4) one can rewrite field equation (3) as 

-fACzz -v) + f,zC P z = o • (8) 

It follows from the equations (7) and (8) the general solution for the metric 
function / takes the following form 

/ = / [C,z - J V(z)dz^) . (9) 
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After introducing new variable £ )Z and derivative /' 

r f df 

Cz = (,z- / n{z)dz, f = -~- 

one can rewrite the derivatives of the metric function / in the following way 



(10) 



f,p — / C,zpj /,pp / C,zp + / C,zpp J 

f,z = f'[C,zz-V(z)], f,zz = f"(<:,zz-V) 2 + f'(C,zzz-V,z). (11) 

Using the equations (7) and (11) one can get 

/,„„ + f,zz + p-'lp = f" [(% + (C,zz - r,) 2 } , (12) 



f 2 P + s% = r [C.p + (c,zz v? 



(13) 



2 , 2 
^■.P+^.z 



c 2 p 2 [C% + (C,zz - v) 2 ] ■ (14) 

Finally, after inserting the derived expressions into the equation (2) we obtain 
equation 

ff n_ f 2 + C 2 f A 







(15) 



which looks like an extension of Papapetrou class of solutions of gravitational field 
equations: 

-l 



f(p, z) = < a ■ cosh 



rj(z)dz 



— (3 sinh 



rj(z)dz 



(16) 



where the functions rj{z) and ((p, z) satisfy to the equation (7) and C 2 — a 2 — (i 2 . 
Here parameters a and (i have the same meaning as in the Papapetrou solutions. As 
a particular case, in the limit r\(z) = 0, the solution (16) reduces to the Papapetrou 



one 



4,5 



3. Asymptotic properties 

As a rule, one of the main difficulties, after finding an exact solution of Einstein 
equations is a physical interpretation of it. In the case of the Papapetrou solutions in 
the Newtonian limit that is at large arguments r = \J p 2 + z 2 — > oo, the asymptotics 
looks like 



f(p,z) 



1+ fc +0 (r- 2 ) 



(17) 



and does not contain a term being proportional tor 1 . From comparison of expres- 
sion (17) with the metric given in 7 

2M 



ds' 



1 



+ 0(r- 2 ) 



dt z - 



S^x k 

4£ijfc— t- + 0(r- 3 ) 



dt ■ dx' 1 - 
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{1 + ^ )6jk+0 (r- 1 ) 



dx j dx k , (18) 



one can conclude that for any given harmonic function £, the metric (16) with i](z) — 
does not contain solution with nonzero mass of the rotating axially symmetric 
source 5 . Here M is the mass of the source and the three dimensional spatial vector 
S J is the total angular momentum. 

We expect that in general case when rj ^ it would be possible to construct a 
solution, which gives flat asymptotics at infinity for massive sources. As an example, 
one of the possible solutions could be obtained in the following way. Assume that 
the function £(p, z) — /i(p)/ 2 (.z) could be taken as a product of two functions /i 
and / 2 . Then from equation (7) one can get 

h,pph + flfazz + P- X h,pf2 = V{Z) • (19) 

Taking solution for the function / 2 in the form / 2 = / 2 , zz = vi 2 ), with rj(z) = 
Co(e z — e~ z ), we have for the function /i the following equation 

OhPP + P' 1 -9h P + 9i -0 . (20) 

Here function f\ = gi + 1 and Co is constant. The last equation is the Bessel one 8 
and /i = 1 + Jo(p), where Jo(p) is zero rank Bessel function. Finally 

t(p,z)=C [l + Jo(p)](e z -e- z ) , (21) 

and according to (10) it gives an expression for the function (, z 

C* = C Jo(p)(e z + e- z ) , (22) 

which takes the simple form C,z|z=o = 2Co</o(p) at the plane z = 0. For the large 
arguments p — > oo, zero rank Bessel function has the following asymptotic behavior 

lim J (p) = J— cos(p - ^) (23) 
p^oo y np 4 

and consequently 

/ 2 7T 1 7T 

C, z |z=o = 2Coa/ — cos(p- -) = — cos(p- -) = x , (24) 
V ^P 4 V P 4 

where we put, for the simplification C = \\f\- Using asymptotic properties of 
hyperbolic functions one could write 

cosh^) = 1 + ip- 1 cos 2 (p - J) + ^p- 2 cos 4 (p - J) + 0(p- 2 ) , (25) 
sinh^) = p-i cos(p - J) + cos 3 (p - J) + 0(p- 2 ) . (26) 
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Substituting obtained asymptotic expressions (25) and (26) in solution (16), we 
can write the asymptotics of function f(p, z), at the plane z — when p — > oo, in 
the following form 



lim f lz=0 = a 1 



1 - \p~ X - \p~' sm(2p) - ^p- 2 cos 4 (p + 



j3 _ i , 7r (3-3 q. n 
+ -p ^ cos(p- -) + — p 2 cos (p--) 
a 4 bot 4 



+ 0( P - 2 ). (27) 



We can see from the asymptotics (27) that in contrast to the Papapetrou one 
(17) there is a term being proportional to p~ x and at the same time we have 
asymptotically flat solution at z — plane. 

4. General solution 

The general solution for function ((p, z) can be derived by the simple integration 
of right hand side of the Poisson like equation (7): 

C(P,*) = -j- r ; ll{Z ' )dZ ' ■ (28) 

Thus in the general case the solution (16) for function / contains function ((p, z) 
which satisfy to the integral (28) and as result function rj(z) generates new solutions 
for the function f(p, z). At the end we would like to underline that for this class of 
functions £, satisfying to the Poisson like equation (7), this is most general solution. 
Using nonvanishing components of Einstein vacuum equations 

Rn = R12 = R22 = (29) 

and following to the book 5 we can easily obtain the connection of function p(p, z) 
with functions ((p, z), u)(p, z) = Cp(^ p and f(p, z): 

^ = -J + JfM & Tp ( < ^ ' (30) 



M,2 — — '-^f + ~T%f,pf,z —U.pUJ.z , (31) 



fjZ P_ r r f_ 

— + -ph P hz - — 

and therefore the metric, presented by expression (1) is defined completely 

One can see form equation (16) that the general solution is invariant under the 
transformation 

C^( + h(z) 

in the functional space of function (, where function h(z) satisfies to the condition 

h(z), zz = rj(z). 

This is gauge condition directly connected with the source of gravitational field de- 
fined by the metric function r](z). From our point of view this statement is analogous 
to the Lorentz gauge condition for the Maxwell equations in electrodynamics. 
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5. Spacetime of line gravitomagnetic monopole 

It is obviously, that any solution of Einstein equations becomes important if it 
has any physical application and identified with the source of gravitational field. 
Here we demonstrate, that the obtained extension of Papapetrou solution contains, 
for example, solution for the line gravitomagnetic monopole and can be useful in 
the interpretation of the connection of the metric source with the gravitomagnetic 
monopole momentum L. 

If we take z) in the form 

j) + j F(z)dz + Vo , (32) 

where m, po, rj are some constants and F(z) — J rj{z)dz, then as it follows from 
the expression (32), so defined function £(p, z) satisfies to the equation (7) for an 
arbitrary function rj{z). Choosing in our solution (16) (i = and a = L/2m, for 
the simplicity, we can obtain 

f(p,z) = ^ — -A — v- ( 33 ) 

L cosh ( 2m In ^ J 

The function (33) coincides with the metric function goo(p) from the paper 9 and 
therefore presents the space time around this source, reproducing cylindrical ana- 
logue of NUT space, named for Newman, Tamburino and Unti 10 . 

From this simple example one can see using (32) that the consideration of the 
line gravitomagnetic monopole metric in the Weyl - Lewis -Papapetrou form leads 
to 

ui = apC tP = 2ma ■ z = L ■ z , (34) 

from where one can conclude, that gravitomagnetic monopole momentum L = 2ma 
in accordance with definition of parameter a. Further investigation of the effect of 
gravitational field of a line gravitomagnetic monopole on test electromagnetic fields 
has been considered in our accompanying paper 11 . 



6. Conclusion 

We conclude, that obtained solution (16) is an extension of the well known Papa- 
petrou class of solutions, which can be applied to the physical systems, presenting 
rotating bounded nonzero masses. We have constructed the solution for the partic- 
ular function r](z), which in some sense, demonstrates the existence of p^ 1 term in 
asymptotics. In addition using properties of Poisson like equation for £ function, 
we have found the general solution (28) for any given function rj(z), which gen- 
erates various metric functions f(p,z). The connection of above obtained solution 
with metric of a line gravitomagnetic monopole was considered as a particular case. 
This is preliminary result and the physical interpretation of the obtained solution 
is outside the framework of this paper. More detailed investigation and application 
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to gravitational sources is the scope of our future research. In this respect as the 
next step we want to investigate the regularity conditions at the axis p = in 
dependence from the choice of function r)(z). 

Acknowledgments 

Financial support for this work is partially provided by the UzFFR (project 01-06) 
and projects F2.1.09 and F2.2.06 of the UzCST. A.K. and B.A. thank TWAS for the 
support under the associate program and IUCAA for warm hospitality during their 
stay there. Authors are grateful to anonymous referee for bringing into attention 
gauge freedom condition, M. Sami for fruitful discussions and Mofiz Uddin Ahmed 
for help and cooperation. A.K. also thanks Yuri Shtanov for the reference to the 
paper 6 and useful comments. 

References 

1. S. Hawking and R. Penrose, The Nature of Space and Time (Princeton University 
Press, Princeton, 1995). 

2. N. Dadhich, arXiv:gr-qc/0311028. 

3. D. Kramer, H. Stephani, E. Herlt and M. MacCallum, Exact Solutions of Einstein's 
Field Equations, ed. E. Schmutzer (Cambridge University Press, Cambridge, 1980). 

4. A. Papapetrou, Ann. Phys., Lpz. 12, 309 (1953). 

5. J.N. Islam, Rotating Fields in General Relativity (Cambridge University Press, Cam- 
bridge, 1985). 

6. P.I. Fomin and V.V. Kuzmichev, Phys. Rev. D. 49, 1854, (1993). 

7. CM. Misner, K.S. Thorne and J. A. Wheeler, Gravitation (W.H. Freeman and Com- 
pany, New York, 1973). 

8. P.M. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw - Hill, New 
York, 1953). 

9. M. Nouri-Zonoz,Class. Quantum Grav., 14, 3123 (1997). 

10. E.T. Newman, L.Tamburino and T.Unti, J. Math. Phys., 4, 915 (1963). 

11. B.J. Ahmedov, A.V. Khugaev and N.I. Rakhmatov, Int. J. Mod. Phys. D 14, 687 
(2005). 



